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Abstract 

This paper presents a conipeUing argument for the physical light speed 
in a homogeneous and isotropic universe to vary with the cosmic time 
coordinate t of Priedman-Lemaitre-Robertson- Walker (FLRW) theory. 
Although this seems to be at variance with the general theory of rel- 
ativity (GR) that assumes a constant light speed, we will show there 
is a reasonable way to accommodate this variable light speed in the 
field equation of GR as applied to FLRW. This variable light speed 
c(t) shows up in a radial transform from a comoving point of FLRW 
to a point on a stationary axis that has a velocity that varies with 
both time and distance due to the time- varying physical differential of 
FLRW. The derivation of c(t) uses the following features: (1) As de- 
rived from the space homogeneity and isotropy of the universe, FLRW 
may have a variable light speed c(t). (2) A generalized time makes 
FLRW metric look the same as constant light speed, so its derivative 
is c(t)dt. This generalized time and its transform back to physical 
time t changes the physics of FLRW to a variable light speed. This 
generalized time when used as the time coordinate of GR also changes 
the physics of GR. (3) The radial covariant transformed coordinates 
describe the same physics as FLRW for radial motion. (4) The coor- 
dinates on accelerated rigid rods of limited length must all measure 
the same c(t), and the same as c(t) of FLRW, a basic assumption of 
relativity theory since all are physical coordinates. Since we assume 
homogeneity of space, this derived c(t) is the physical light speed on all 
points of the FLRW universe. This has impact on the interpretation 
of all astronomical observations of distant phenomena. 
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1 Introduction 

This is a significantly different approach than other attempts in the liter- 
ature to investigate a variable speed of light. Those mostly tried to find a 
new cosmology to provide alternatives to inflation in the "standard" FLRW 
cosmology to resolve horizon and flatness problems [3] [4]. 

I did not start with an intent to find a variable speed of light. I was 
interested in finding a connection between the stationary rulers and clocks 
in our galaxy with the FLRW co-moving coordinates of a homogeneous and 
isotropic universe. As standard texts on general relativity make clear, it 
is not possible to find a coherent theory of gravitation using three rigid 
orthogonal axes in an accelerating environment. For this we need the curve- 
linear coordinates of general relativity. But it did seem possible to find a 
radial transform from the comoving FLRW coordinates to the local physical 
coordinates. This led to the discovery of the variable physical light speed. 

The FLRW metric[l] [2], which can be derived from the assumption of 
homogeneity and isotropy, allows for a variable physical light speed c{t) that 
depends on the cosmic time t. To display this we can use a generalized time 
t, where dt = c{t)dt, to make the FLRW metric have the same coefficients 
as for constant light speed. Of course, using t and its transform back to the 
physical time t make a real change in the physics because it implies that the 
light speed varies with physical time throughout the universe. 
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When one set of physical coordinates transform to another, the physics 
doesn't change, but the transform can be constrained. Thus, when the ordi- 
nary coordinates on a rigid frame are transformed to another that is moving, 
and insist that both measure the same constant light speed, this constraint 
gives the usual Minkowski metric and the Lorentz transform. When these 
are allowed to have the same variable light speed, this gives the generalized 
Minkowsli metric (M) and generalized Lorentz transform {dL). The FLRW 
comoving coordinates i, x have a metric like M for radial world lines, for 
which dL covariantly transforms FLRW physical radial differentials to the 
colocated differentials of radial inertial rods of limited length. Because of the 
time variation of the FLRW physical radial differential a{t)dx, the velocity 
between a galactic point and a point on the stationary axis is a function 
of both time and distance. This dual dependence places an additional con- 
straint on the transform that determines the function c(t). 

This is shown by expanding the physical time and distance along a sta- 
tionary rod in a power series of the FLRW distance x from the origin using 
dL from the FLRW coordinates. Although not a priori obvious, in the limit 
of zero distance while keeping the lowest significant power of x, these con- 
straints lead to c{t) oc y^da/dt. 

The reasonableness of this limiting process is confirmed by full radial 
covariant transformed coordinates from FLRW, good for all distances, whose 
differentials have a M metric close to the origin. These all have the same 
variable light speed at the origin shown by the power series expansion. 

For a homogeneous universe, since the origin can be placed on any galac- 
tic point, this means that this variable physical light speed enters all our 
physical laws throughout the universe. 

We can accommodate the variable light speed c(t) in the field equation 
of general relativity for FLRW by interpreting the time variable to be t in 
the reduced curvature tensor and keeping the gravitational and rest mass 
energies constant in its proportion to the stress-energy tensor. This is a 
real change in the physics, and not just a coordinate transform, because 
t has changed the physics of the FLRW line element. This enables us to 
calculate a{t), c{t), and trajectories in the time and distance of stationary 
coordinates. 

We can express the gravitational field in the transformed coordinates. In 
the region near the origin for a flat universe this held increases linearly with 
distance just like the Newtonian field for a spherical distribution of uniform 
mass density. 

A surprise bonus from this endeavor is that one of the radial covariant 
transforms has a physical distance to all parts of the universe. Even though 
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three rigid accelerated axes are inadequate to describe three dimensional 
motion, it is apparently possible to find one rigid axis to measure radial dis- 
tance, at least for a homogeneous FLRW universe, although the transformed 
time on this axis becomes non-physical at large distances. 

2 Assumptions 

Because a variable light speed is such a major departure from the accepted 
paradigm that light speed is constant for all times throughout the universe, 
I will first list the assumptions necessary for the derivation of c{t). These 
are the usual assumptions for an FLRW universe adapted to the possibility 
of a variable light speed. The assumptions in quotation marks are enumer- 
ated with italic Roman numerals that will be used throughout the paper 
when the assumption is being used. The remarks not in quotation marks 
are descriptions, definitions or derivations. Assumptions i — v arc needed to 
find c{t) and vi — viii are used to apply c{t) to physical situations. 

i. "The universe is isotropic and homogeneous in space." 

The generalized FLRW metric (eq 6) is derivable from this Cosmological 
Principle (Sect 3) with ds = c{t)dT as the generalized invariant line element. 
The adjective "generalized" will be used to indicate in all equations that c{t) 
as found in the FLRW metric has replaced a constant cq; in particular that 
the generalized time t, for which d,t = c{t)d,t has replaced the cosmic time 
Cot [5]. (t has the dimensions of length). It should be emphasized that this 
line element comes from the assumed symmetry and not as a solution to the 
field equation of general relativity. As such, the use of t in the line clement 
and its transform to physical t will create different physics in FLRW char- 
acterized by a variable light speed. 

a. "The usual interpretation of the FLRW coordinates applies to the gen- 
eralized FLRW coordinates." 

The line element ds for the FLRW metric has time and radial distance 
differentials of dt, a{t)dx, that represent physical quantities, measurable by 
clock and ruler differentials or their technological equivalents, so that the 
physical light speed is c(i) = a{dx/dt)s, and the physical radial velocity 
of a moving object, labeled R, located ai t,x '^s a{dx/dt)ji. A galactic 
point is denoted by %, so x is a co-moving radial coordinate with which a 
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galactic point stays constant, t is the time on a clock co-moving with each 
galactic point, including the origin at x = 0. Each galactic point is on a fi- 
nite inertial frame of limited extent, which moves with it. For radial motion 
the FLRW metric becomes Minkowski- like with a variable radial differential. 

Definitions: We will use "covariant" to describe coordinate systems x'^ 
that are covariant transforms from the FLRW coordinates. They thus have 
an invariant line element ds equal to the FLRW line element and a con- 
travariant velocity vector dx^^ /ds (dimensionless) whose coordinate compo- 
nents transform like the transform of coordinate derivatives. We will use 
the word "full" to describe a coordinate transform for the full range of co- 
ordinates, and "approximate" to describe a transform for a limited range of 
coordinates. 

We will consider a radial world line measured with an isotropic {i) ra- 
dial covariant transform (Sect 6) from FLRW {ii) whose time and radial 
coordinates are T{t^ x) ^-iid ^{ii x) ^iid whose space origin is attached to 
the same galactic point as x = 0, at which point there is no motion be- 
tween them, which we will therefore call a stationary frame. We will use the 
word "colocatcd" when the coordinates on two different frames are at the 
same space-time point. We define a generalized Hubble ratio as H{t) = d/a, 
where the dot is the t derivative. 

in "Transforms of the colocated differentials of physical coordinates be- 
tween accelerated inertial frames are generalized Lorentz transforms {dL ) 
with the same physical c{t) as FLRW (Sect 3) in the limit of small frames. 

It is well known and experimentally verified for physical coordinates on 
earth in the present era, that the transforms of special relativity apply lo- 
cally to inertial frames in an accelerating field. It seems reasonable that 
this could be generalized for a different cosmic time which might have a 
different physical light speed. The differentials of physical coordinates on 
these inertial frames will therefor be represented by a generalized Minkowski 
metric M (Sec 4). Both dL and M can be derived from the requirement 
that all moving observers measure the same light speed c{t). Although the 
transform differentials are not integrablc to full Lorentz transforms between 
inertial systems in free fall, they will be integrable for small ranges of time 
and distance. We will investigate how large a range will reasonably repre- 
sent physical coordinates (Sec 6). 

(iv) "Close to the origin, the stationary frame is inertial and has physical 
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coordinates whose differentials have an SR M metric and for radial world 
lines are dL transformable from the physical differentials of the FLRW met- 
ric. " 

Since the origin is located on an inertial galactic point, the frame at- 
tached to it must be inertial close to it and have differentials of physical 
coordinates {in). 

(v) "The r, i? of a radial covariant transform from t, x whose differentials 
have an M metric close to the origin represent physical coordinates (iv) in 
the same range that the transformed metric is M ." 

Although it is conceivable that there are transformed coordinates from 
t, X with M metric that are not physical, the fact that all such found systems 
have exactly the same c(t) close to the origin reinforces this assumption of 
their physical uniqueness. It is also reinforced by the Bernal criteria (vi) 
which show that the differentials of the dL transformed coordinates use the 
same clocks and rulers as do the physical differentials of FLRW (ii) close to 
the origin (see section 6). 

With these assumptions and definitions, we will show that the variable 
light speed c{t) is proportional to d = aH (or equivalently to ^Jda|dt) by 
two different procedures: 

1) Integrate dL transformed physical differentials dT, dR for small x {iv). 
(Sections 4 and 5) 

2) Find full diagonal radial covariant transforms T, R for all t, x that 
have a M metric for small x (v)- (Sections 6 and Appendices A-C) 

Each of these has the same light speed in the limit of X^O. The first 
shows this for any and all covariant transforms that are physical close to 
the origin for an expansion of T to the second power of X: the first order 
vanishing. The second shows this for a large number of full radial covariant 
transforms which have a M metric close to the origin. Thus, the first is a 
completeness proof that if there are such transforms, they must have this 
c{t), and the second is an existence proof that there are such transforms and 
that the expansion of the first is further justified for being carried out to 
second power of %. 
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Additional assumptions are needed to apply this variable physical light 
speed to physical laws. A possible approach is the following: 

(vi) "We assume the Bernal criteria[12] that two observers will be using 
the same units of measure when each measures the other's differential units 
at the same space-time point compared to their own and finds these cross 
measurements to be equal." 

We will find a radial covariant transform (T, R) called physical distance 
coordinates whose difi^erentials dR are all physical by virtue of this assump- 
tion (Section 6 and Appendix B). If wc assume that thereby the integrated 
i? is a rigid physical distance, this R becomes a global reference axis (Section 

6) . 

{vii) "We assume that the time in Einstein's field equations can be in- 
terpreted as t, where dt = c{t)dt, and that the rest mass energy is kept 
constant." 

For FLRW this makes the usual da/dt in the field equation become 
da/dt = a for a variable c{t) (Sects 7 and 8). The solution for t converts the 
physics to a variable light speed. It's transform to the physical t provides 
physical coordinates to describe this physics. This allows us to calculate 
galactic and photon paths for a homogeneous and isotropic universe with a 
variable light speed. In addition, we assume the rest mass m of a particle be 
proportional to 1/c^ to keep the rest mass energy constant. The field equa- 
tion then requires the gravitation "constant" G to be proportional to (Sect 

7) . This keeps constant the static gravitational energy. The geometrization 
of special relativity can be preserved with t^t, and m^m = mc^ [?]. 

Although we prove c{t) oc a by using special relativity applied to the 
FLRW metric, the field equation shows that gravitation caused by the uni- 
verse energy density determines a{t) and thus c(t). We show that for a flat 
homogeneous universe there is a gravitational field that increases with R 
just like the Newtonian field for a spherical distribution of uniform mass 
density (Section 9). 

{yiii) We assume that the fine structure constant and Rydberg frequency 
remain constant. " 

This seems to be required in the present era by current measurements 
with atomic clocks [6] [7]. This is accomplished by assuming the vacuum 
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electric constant eo{t) and the vacuum magnetic constant /xo(t) each to be 
inversely proportional to c{t) (Sect 10). With this assumption, the defini- 
tions of the electro-magnetic fields can be generalized to preserve Maxwell's 
equations and the electro-magnetic field tensor. 

Some suggestions are made in Sect 10 on how c{t) might be observed. 



3 The generalized FLRW metric 



Implicit in i is that the concentrated lumps of matter, like stars and galaxies, 
can be averaged to the extent that the universe matter can be considered 
continuous, and that the surroundings of every point in space can be assumed 
isotropic and the same for every point. 

By embedding a maximally symmetric (i.e., isotropic and homogeneous) 
three dimensional sphere, with space dimensions r, 9, and (j), in a four dimen- 
sion space which includes time t, one can obtain a differential line element 
ds [8, page 403] such that 



ds^ = g{t)df - fit) 



dr"^ 



1 — kr"^ 



+ r'^de'^ + r'^sin'^edf 



where 



r = < 



sinx, A: = 1, 
X, A; = 0, 
sinhx, A; = — 1, 



(1) 



(2) 



A; is a spatial curvature determinant to indicate a closed, flat, or open uni- 
verse, resp., and 

dx^ = dr'^/{l-kr'^). (3) 

Wc let a{f) = ^/Jit) be the cosmic scale factor multiplying the three 
dimensional spatial sphere, so that the differential radial distance is a{t)dx- 

The g{t) has normally been taken as g{t) = (? = constant, so that c is the 
constant physical light speed and t is the physical time on each co-moving 
point of the embedded sphere. In both cases by physical, we mean that their 
value can represent measurements by physical means like standard clocks 
and rulers, or their technological equivalents. In order to accommodate the 
possibility of c{t) being a function of time, we make g{t) = c{t)'^. The 
resulting equation for the differential line element becomes a generalized 
FLRW metric: 

ds^ = c{tfdt'^ - aitfldx^ + r'^diu% (4) 
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where doS^ = dO'^ + sin^0d4>^. For radial world lines this metric becomes 
Minkowski in form with a differential of physical radius of a{t)dx- 

It will be convenient to introduce the time related quantity t [5], which 
we will call a generalized cosmic time, defined by 



where c(f) = c(t), and where the lower limit is arbitrarily chosen as 0; this 
choice does not affect the derivation of c(i). 
The line element then becomes 



It should be emphasized that i itself is a legitimate more general coordi- 
nate that makes the metric of FLRW satisfy the Cosmological Principle; the 
physical time t is a transform from it. t and its transform to t change the 
physics to a variable light speed. This is useful in GR since the coefficient 
of the differential of i in the FLRW metric is constant just like for t with a 
constant light speed. 

Notice that this invariant line element is different from another line cl- 
ement dr = ds/c{t). For a constant light speed, these two elements are 
interchangeable, but for the variable light speed being considered here, they 
are not. r is still the time on clocks in freefall in the universe, like the time 
on cosmic points x, but is not taken as invariant under transforms. 

4 Generalized Lorentz transform from galactic points 
to the stationary inertial frame using the veloc- 
ity V between them 

We will consider only radial world lines with physical coordinates T and R 

on the stationary inertial frame. We would like these to describe the same 
events as t and X) so T = T(i, x) and R = R{t, x) with i2 = at x = 0. So 



where the subscripts indicate partial derivatives with respect to the subscript 
variable. 





ds^ = dp - a^idx^ + r^dw^). 



(6) 



dT = Ttdt + T^dx = \Ttdi + T^dx, 
dR = Rtdt + R^dx = \Rtdi + R^dx, 



(7) 
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We will find T = T{t, x) and R = R{t, x) by integrating the differen- 
tials of the Lorentz transform. We have extensive measurements for phys- 
ical systems on earth that confirm that the Lorentz transform of difi'eren- 
tials between inertial frames is locally accurate and invariant, and keeps 
the Minkowski metric invariant when c{t) is substantially constant, even in 
the presence of a gravitational field. If the physical light speed is going 
to change with time, it seems reasonable that physical systems retain the 
generalized Lorentz transform of difi'erentials dL between inertial frames of 
limited size using the physical light speed c{t) , which would keep the gener- 
alized Minkowski (M) metric invariant on any moving inertial frames (iv) 
for small frames: 

ds'^ = c{tfdT^ - dR^ - R^du'^ , (8) 

This M metric (eq 8) is invariant under an invariant dL between a second 
inertial frame {dT, dR) moving at a radial velocity V{t) relative to a given 
inertial frame {dT',dR'): 

dT = ^{t){dT'-y^dR'), 
dR = -f{t){-V{t)dT' + dR'), 



where j{t) = l/y 1 — V{t)^/c(t)^. We assume the M metric applies to phys- 
ical diflcerential times and distances on inertial frames of limited size any- 
where and anytime {iv). The FLRW metric in eq 4 has a radial Minkowski- 
like metric with dT'^dt and dR'^adx- If a point on the stationary firame 
is moving at a radial velocity V{t, x) when measured with the FLRW coor- 
dinates, the dL transform of dt, adx to dT, dR for a radial path keeps the 
line element ds invariant: 

dT = ^{t,x){dt-^a{t)dx), ^^^^ 
dR = 7(t, x)i-Vit, x)dt + a(t)dx), 

The basic assumption of the FLRW cosmology is that the differentials dt 
and adx of the FLRW metric arc physical. 
If we compare eq 10 with eq 7, we get 

Tt = l, (11) 
= -jaV/c^ = --raV/c, (12) 
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Rt = —fV = -jcV, 



(13) 



Rx = 70, (14) 

where for simplification we have introduced V = V/c. These relations are 
exact for differentials as x^O, and therefor are approximately correct when 
the differentials are integrated for small We can rearrange the two ex- 
pressions for V to give 

With cqs 11 and 14 this gives two relations each for dT and dR in terms of V . 
When we try to integrate these partial differential equations, we integrate 
dT, dR along the R frame but integrate the dt, dx along a radial connection 
between the comoving galactic points x- Because the radial differential 
changes with time because of a{t), V{t,x) changes with time and distance. 
As we integrate dR with the time dt at constant x, cV must vary in such a 
way as to maintain dL . Simultaneously, as we integrate T with x at constant 
t, aV /c must also vary to maintain dL . The combination of these double 
constraints requires c(t) to vary with t in a determined way, at least for the 
short distance from the origin where dL is valid. We will find that if there 
is no acceleration, and V is a function only of x ™ an expanding universe, 
this removes one constraint, and c{t) will be constant (see Appendix A. 5). 



5 Variable light speed c(t) derived from the Lorentz 
transforms 

To obtain T(t, x) and i?(t, x) near the origin, we need to integrate the differ- 
entials dT and dR for small %• We will do this by expanding these physical 
coordinates in a power series in x out to the lowest power that will give a 
non-trivial c{t) in the limit of zero x- We will use the two relations for dR 
to determine the expansion coefficients of R and V, then use the resultant 
expansion of V in the two relations for dT to expand T and determine the 
requirement for c(t). 

Since V will vanish at the origin (see definitions. Sect 2), the constant 
in the power series for V is zero; so let 

-V = wi{t)x + W2{t)x^ + ws{t)x^ + Oix")... (16) 
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where the Wi{t) are unknown functions to be determined. Prom = 07 
(eq 14) we get 

= a(l + + ^V^ + ...) = a(l + ^-wjx' + W + 0(x'))- (17) 

If we integrate eq 17 at constant t, noting that R vanishes at % = (see 
definitions, Sect 1), we obtain 

R = ax+ ^awjx^ + ^awiW2X^ + 0(x^)- (18) 

R{t,x) is the physical differential dR summed over all the galactic points 
up to Xi s-iid is thus the physical distance to x a-t time t. The first term of 
eq 18 is the "proper" distance to which all measurements of distance reduce 
close to the origin [8]. 

Partial differentiation of eq 18 by t at constant x gives 

Rt = cax + Ix'^^^iawj) + ^x" j^{awiW2) + 0(x^)- (19) 

where the dot represents the derivative with respect to i. We can then find 
V from eqs 15, 17, and 19: 

-V^ = 1^ = OX + f{t)x' + 0(X')-, (20) 

where 

By comparison of eq 20 with eq 16, we see that u^i = a, W2 = 0, and 

W3 = fit). 

We will now use this expression for V to find two relations for Tj. The 
first comes from Tt = 'j (eq 11): 

Tt = l+^-V' = l + \a\^ + 0(x')-, (22) 

Even though dT and dt arc both measured on standard physical clocks, we 
note that the galactic clocks t measured at constant x I'un slower than the 
stationary clocks T as they move away from the origin (dilation). When 
measured at constant R, the stationary clocks run slower, tT = l + , in 
accordance with the Lorentz transform. Of course the distance contraction 
is also consistent with Lorentz, Rx/a = aXR = 1 + V^^/2 (eq 17). 
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We can find an expression for T^, using eqs 15, 22 and 20: 

= --TtV =-[! + \i?x^ + 0(x')-][ax + /Wx' + 0(x')-], (23) 
and multiplying the brackets gives 

= "[«X + + /Wx' + 0(X')-]- (24) 

By integration with x constant t with T = t at x = (ii) we find 

r = i+2fx' + 0(x')- (25) 

We have thus T(i, x) at i?(i, x) requiring a synchronization of the physi- 
cal stationary clocks that changes with time either at constant x or constant 
R. This should not be surprising in special relativity for a free falling rigid 
rod in a time varying gravitational field that is inertial at the origin that 
becomes non-inertial with distance along the rod. 

If we partially differentiate eq 25 by we get a second expression for Tj: 

The equivalence of the two expressions for Tt (eqs 22 and 26) out to the 
lowest non-trivial term in % requires that there be a variable c(t) given by 
the diflFerential equation: 

Mathematically, when we regard c(t) as a variable to be determined by 
the limiting process of x ^ 0, we must keep the term in since it is 
the lowest term that determines c(t), which we have therefor called non- 
trivial. Remember that the c(t) in these equations is the physical light 
speed assumed for the FLRW metric (ii) and for the M metric (iv). 

To get an explicit expression for c(t), multiply eq 27 by a, change the 
variable dt to da/cd to yield 

^ = ^ci(^). (28) 
a aa c 



One can see that c oc d is a solution, so 

c{t) a{t) 



Co d(io) 



OiE, (29) 
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where a is the normahzed scale factor 



a = a/ao, (30) 
and E is the normahzed Hubble ratio 

The subscript denotes the value at t = to, the present time. We can take 
Co to be unity, so that c(t) would be measured in units of cq, but for most 
equations in this paper I will retain cq for clarity. The field equation (sect 
7) will enable us to evaluate a and E and thus c{t). 



6 Radial covariant transforms that become physi- 
cal close to the origin 

6.1 Procedure for finding radial covariant transforms using 
the velocity V 

We would now like to find radial covariant transforms that will hold for all 
values of the FLRW coordinates and reduce to the physical coordinates for 
small distances from the origin. The most general line element for a time 
dependent spherically symmetric (i.e., isotropic) line element (Weinberg, 
p335) is 

ds'^ = c^A^dT'^ - B^dlW- - IcCdTdR - F'^{dd^ + sin^ed(^^) (32) 

where B, C, and F are implicit function of T and R, but explicit functions 
of t and X- We are using the same notation for time and distance as we did 
for the physical coordinates, but understand that they may be physical only 
for small distances from the origin. We have included the physical light 
speed c{t) in the definition of the coefficients of dT. 

We will look for transformed coordinates which have their origins on the 
same galactic point as x = 0, so i? = when x = 0, where there will be no 
motion between them, and where T is t, since the time on clocks attached to 
every galactic point is t (ii), including the origin. We will make F = ar to 
correspond to the FLRW metric, but will find only radial transforms where 
the angular differentials are zero. Of course, full four dimensional transforms 
to time and three rigid axes have not been found, nor are they required to 
determine c{t). They don't even have to be covariant. They have only to 
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meet the requirement of becoming dL close to the origin. Radial covariant 
transforms are convenient to do this. Then R and T will be functions of 
only t and x- T = T{t, x) and R = R(t, x), and we will still have eq 7. Let 
us consider a radial point at R in the transformed system. When measured 
from the FLRW system, it will be moving at a velocity (ii) given by 

This velocity will be the key variable that will enable us to obtain radial 
covariant transforms of the full coordinates. To get the time component of 
the contravariant velocity vector = di/ds, we divide eq 6 by dP with 
duj = to obtain 



(|)^.i_„,)^(|)^.a_f^).J, (34) 



We can make this substitution of dx/dt = V/a from eq 33 because we are 
calculating the transform for the point R. To get the spatial component, we 
use the chain rule applied to eqs 33 and 34: 

dx dxdi V 

— = —— = —7. (35) 
as dt ds a 

The spatial components of the contravariant velocity vector of a point 
in the transformed system are zero because we will assume the point is not 
moving in those coordinates. This assumption means that a test particle 
attached to the radial coordinate will feel a force caused by the gravita- 
tional field, but will be constrained not to move relative to the coordinate. 
Alternatively, a colocated free particle at rest relative to the radial point 
will be accelerated, but will thereafter not stay colocated. This means that 
the radial point is not inertial. However, it also means that the radius will 
be an axis of simultaneity since there is no motion between its points. (By 
simultaneity I mean that any two events that are simultaneous at one point 
on R are seen as simultaneous at all points via light signals). 

The time component of the transformed system is dT/ds = 1/cA. This 
makes the vector 11^^ = in the transformed coordinates 

C/'^ = (0,0,0,^) (36) 



and in the FLRW coordinates 



= (^,0,0,7). (37) 
a 
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Since it is contravariant, its components will transform the same as dT, dR 
in eq 7: 



cA 



= lR,j+lR^jV. 
Manipulating the second line of eq 38 gives 

aRt 



V 



cR^ 



If we invert eq 7, we get 

di= ^{\R^dT -T^dR), 
dx= U-c^tdT + TtdR), 

where 



D = -TtR^ - -RtT^ = -TtR^{l + F^), 
ccc alt 



(38) 



(39) 



(40) 



(41) 



using eq 39. Because of the invariance of ds^, we can enter dt and dx of eq 
40 into the FLRW metric (eq 4 ) to find the following expressions for the 
coefficients of eq 32: 



A' 



and 



C 



1 



a 

TtRy 



(1+Fgl)2 



If we put ds = in eq 32, we obtain a coordinate velocity of light Vp. 

'dR] 



dT 



(42) 

(43) 

(44) 
(45) 



We need to remember that the c{t) in these equations is the physical light 
speed assumed for the FLRW metric {ii). The equations for A, 5, and Vp 
simplify for a diagonal metric (C = 0): 



A 



1_ 
Tt 



tr 
7 



(46) 
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where we have used eq 40 with C = to obtain the inverse partials. 

This metric becomes M when A ^ 1,B ^ 1,C— ^0 and ar — ^ R, 

and we get the relations in eqs 11-14. Even when these conditions are 
not met, we can say something about the physicality of the coordinates 
with the use of criteria (vi) developed by Bernal et al[12]. They developed 
a theory of fundamental units based on the postulate that two observers 
will be using the same units of measure when each measures the other's 
differential units at the same space-time point compared to their own and 
finds these cross measurements to be equal. Thus, even ii B 1, dT will be 
physical if C = and A = 1 because then Tt = tr = J so that dT uses the 
same measure of time as dt, which FLRW assumes are physical. Similarly 
even if A / 1, dR will be physical if C = and B = 1 because then 
Rx/o- = o^xr = 7 and dR uses the same measure of distance as adx, which 
FLRW assumes are physical. If we assume R is physical when physical 
dR is integrated indefinitely out to the visible horizon, since there is no 
movement between any two points on i?, it becomes the radius of a rigid 
global reference frame of simultaneity. Although it is not possible to find 
physical global coordinates in four dimensional frames in a gravitational field 
([9], it apparently is possible to find a global physical radius. 

At this point we would like to examine quantitatively how far from the 
M metric our transformed metric is allowed to be in order for its coordi- 
nates to reasonably represent physical measurements. We can consider the 
coefficients A, B, and C one at a time departing from their value in the M 
metric. For example, let us consider the physical distance case B = 1,C = 
and examine the possible departure of the time rate in the transform from 
that physically measured. Then, from eqs 46: Tt = j/A,tT = jA. Thus, 
1 — A represents a fractional increase from 7 in the transformed time rate 
Tt, and thus the fractional increase from physical {v) of an inertial rod at 
that point. We can make a contour of constant A on our world map to give 
a limit for a desired physicality of the transform (see Figs 1-4, where the 
heavy dots represent 1 — A = 5%) . 
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6.2 Diagonal radial covariantly transformed coordinates with 
physical c(t) close to the origin 

We show in Appendices A, B, and C that there exist an infinite number 
of radial covariantly transformed coordinate systems which satisfy the M 
requirements close to the origin. Appendix A derives diagonal transforms 
(C = 0) using physical time {A = 1) for all distances R, which all show that 
the light speed becomes c{t) oc a when the transforms becomes Lorentz. 
Appendix B shows the transforms for physical distance {B = 1,C = 0) for 
all times T. Appendix C describes similarity solutions of these equations 
for a flat universe (O = 1). Because the equations in A and B are different 
from each other, they show, as we would expect, that it is not possible to 
have diagonal transforms with physical R and physical T simultaneously for 
all values of t,x (except for an empty universe). Wc can thus picture the 
universe with stationary covariant coordinates extending out through the 
universe that is physical with a dL transform from other co-located galactic 
points for a small distance from the origin. 

At all distances for A = 1, the covariant time T can be measured on 
continuously synchronized stationary physical clocks, but the covariant dis- 
tance R cannot be measured on physical rulers for all distances. For B = 1, 
the covariant distance R can be measured by physical rulers on a station- 
ary frame for all distances, but the covariant time T cannot be measured 
by physical clocks (except for small V). We can calculate an acceleration 
(Appendix D) for a fiat universe that is zero at the origin, and increases 
with distance; the physical distance R acts like you might expect for a rigid 
ruler on whom the surrounding masses balance their gravitational force to 
zero at the origin, but develop an inward pull as the distance increases. 

Because there is a special interest in having a physical description for 
distance in the universe, we display the physical distance transforms. The 
physical distance results for flat space (Appendix C) are shown in Figs. 1 
and 2. Here we have used the field equations with the generalized time (Sect 
7) to derive the equations for a{t) = aoit/toY^'^ and c{t) = co{to/t)^^^- 

Fig 1 plots distance R against the time at the origin (cosmic time t) for 
galaxies (constant x) ^ov incoming light reaching the origin at t/tQ = 1. 
The galactic paths are labeled with their red shift z, determined by the 
time t of the intersection of the photon path with the galactic path z = 
— 1 -|- c/a = — 1 + {to/t)^^^, assuming the frequency of the emitted light 
does not change with c{t). Notice that light comes monotonically towards 
the origin from all galactic points. This photon path has a slope of cq = 1 
close to the origin where the distance R and time t are both physical, but 
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decreases as the distance increases and the time decreases, different from 
c{t). 

Although the distance uses physical rulers, the coordinate system as a 
whole may not be physical for times shorter than some limit. A reasonable 
limit might be ^ = 0.95, R = 2.3{t/to)^^^, shown by the heavy dotted line 
in the figures. Together with B = 1,C = for these physical distance coor- 
dinates, this shows that the assumption that T and R represent a physical 
stationary coordinate system inside this limit is very good, with distances 
accurately represented, and time rates Tt within 5% of physical measure- 
ments on adjacent inertial rods. 

Fig 2 plots these distances vs the transform time T at R. At the emission 
of the photons, T is finite (even for t = 0), presumably the transformed 
time it takes for the galactic point to get out to the point of emission. At 
T/tQ = 1 the slope of the light path is cq = 1, and at the physicality limit 
T/tQ = 0.40 the slope is only 5% less than c(t) = 1.50. At the intersection of 
this physicality limit with the photon path that arrives at the origin at to, 
the time t/to = 0.2 and the red shift z = 2.4. Thus, if we have a flat universe 
with 17 = 1, the last 80% of the universe history out to a z of 2.4 can be 
treated with physical coordinates T and R. This z is as large as any of the 
Supernova la whose measurements have suggested an accelerating universe. 
Even out to z = 7, ^ - 1 = 0.2, = All, B = 1,C = 0, Vp/c = A = .8, 
c{t) = 2, and t/to = -06, so the approximation to M is not too bad for a 
universe 6% of its present age, with a light speed at twice its present value. 
The locus of points where t = 0, labeled horizon, is interesting, but does not 
represent physical coordinates (see Sect 9) . 

Actually, the coordinates of the transform are closer to M than keeping 
only the lowest order of the expansion. Thus, for the plotted physicality 
limit, z = 2.4, keeping the lowest order of the expansion of Tt gives A = 
'y/Tf = 7/(1-1- F^/2) = 1.14, because the lowest order of the expansion is 
less than 7; whereas for the transform, A = 0.95, making the transformed 
Tt a slight overcorrection, but much closer to 7. 

I have also included three additional figures, also using the generalized 
field equation of section 7 and 8. Fig 3 is for a density of = 1/2 (Appendix 
B), which has paths intermediate between Q = 1 and $1 = 0. Fig 4 shows 
the effect of dark energy (Appendix B for JIa = 3/4), where all the curves 
tend to have infiection points when the dark energy becomes dominant. 
The empty universe (fi = in Appendix A. 4) shown in Fig 5 is physical for 
all space-time, undistorted by gravitational curvature; galactic points and 
light travel in straight lines. Figs 1-2 are from the numerically integrated 
similarity solution (App C.l), Figs 3-4 are from the numerically integrated 
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initial value solution (App B), and Fig 5 is an analytic function solution 
(A. 4. These illustrate complete coverage of < O < 1. 

7 Revised Equivalency Principle and Einstein field 
equation 

A number of different ways to accommodate a variable light speed have 
been addressed in the literature [3] [4]. In order to make calculations of the 
present determination of c(t), we are here adopting a simple accommodation. 
If this should prove to be wrong, it still should not invalidate the derivation 
of c{t) oc a. To accommodate c{t), it would seem simplest to revise the 
Equivalency Principle as stated by Misner et al[9, page 386] to read "to 
any and every local [generalized] Lorcntz frame, anywhere and any time in 
the universe, all the (nongravitational) laws of physics must take on their 
familiar [generalized] special relativistic forms {in)", viz. dL (eq 9) and the 
M metric (eq 8). 

The use of di[5] in the field equation (vii) seems to be consistent with 
this: 

87rG' 

G^^ + A<7^, = ^ V, (49) 

where G^;/ is the Einstein contracted curvature tensor determined from the 
coefficients Qui, multiplying the differentials in the metric of eq 6, A is the 
cosmological "constant" possibly representing some kind of vacuum energy 
density, G is the gravitational "constant" and T^^ is the energy-momentum 
tensor, c has been inserted in the ratio between them in order to balancee 
dimensions. 

It will be convenient to use the variables t, r, 6, <p in the generalized FLRW 
metric (eq 6): 

ds"^ = dp - a2(- — + r^de^ + r'^cos^edcp'^). (50) 

so the gfj,i, of eq 49 are the coefficients of eq 50, exactly the same as the usual 
FLRW metric with constant c = 1 in the variables t, r, 6, (j). 

It should be emphasized that t is the solution of the field equation repre- 
senting the physics of a variable light speed. Physical time t can be acquired 
by transforming from t (eq 5), preserving this physics. 

We will make the usual assumption that the universe is an ideal fluid 
with an energy density of pc^ and pressure p, so that we can write the two 
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significant field equations [9, page 729] for a{t) as 



3d^ 3k , SttG 

-^ + ^-A = ^p, 51) 

and 

a A: . SttG , , 

+2- + ^ + ^-A = ^p, 52 

where the dots represent derivatives with respect to t. 

Now, the first equation can be multiphed by a^/3, differentiated, and 
subtracted from ha} times the second to give 

d Gp 3 3G . 2 . 

^(^«) = -^-^'-^A, (53) 

where we have allowed for the possibility that G, p, p, and A may be func- 
tions of c(i). For small p and A, 

Gpa^/c^ = constant (54) 

If the energy density consists of n particles per unit volume of mass m, so 
p = nm, then the conservation of particles requires na^ be constant. This 
makes 

Gm 

— = constant. (55) 

Even when p and A may not be negligible, it is reasonable to assume that 
for matter density eq 55 will remain valid. The relation of constants in 
eq 55 also preserves the Schwarzschild metric as a generalized solution of 
the field equation, since it requires for its derivation that the Schwarzschild 
radius 2Gm/c^ be independent of time. Consistent with the conservation of 
the stress energy tensor, we assume m oc and by eq 55, G oc c^. This 
keeps constant both the rest mass energy and the static gravitational energy. 
Some implications of this on particle kinematics and on electromagnetism 
are treated in Ref [?] . 

8 Calculation of a(t) and c(t) 

Following Peebles [10, p312], we define 

_ SttGo , . 
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and 

^r^^, (57) 

and 

i^A = ^. (58) 

For very small a there will also be radiation energy density p = (cq /cj'^gaoT'^, 
where a^T^ is the black body energy density of the various relativistic par- 
ticles represented by g. If T oc 1/a with no c{t) dependence, we can define 
a radiative term cofib/ca^, where 

n, = gao{Ta)'^^. (59) 

The normalized Hubble ratio £^ in eq 31 is determined by eq 51: 



Id CQ^lb ^ ^ 



= £;=W^+ +^ + 0a. (60) 
which allows us to evaluate c(t)/co = aE. The Jls are defined so that 

nb + n + nr + nA = i. (6i) 

At t = to: a = 1, E = 1, and c/cq = 1. 
The cosmic time t becomes 

For J7 = 1 (J^b = = = 0): 

coi^oio = 1/2, (63) 
a = (t/to)^/' 

For other densities with = = 0, 

co^oi= ^^4^b-ln(l + y)], (65) 



where 



y = (66) 
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There is no periodicity of a with t for O > 1, . The higher density de- 
creases the time to^ln^/coHoQ,, and the universe scale factor a contin- 
ues to expand, asymptotically approaching a maximum at $7/(17 — 1). As 
ri-H'O, a-^i, c-^co, the universe becomes Minkowski (see Appendix A. 5). 

The radiative term becomes ft^^^/a^^^^ for a « O^/Oq- ^6 is of the 
order of 10"^ [10, pl64] and only dominates the expansion rate when 1 + z = 
c/coa = E » il^/O^f, which is of the order of 10"^^. The inverse c{t) 
dependence of the radiative terms reduces its influence on a and extends 
the mass dominant era to much shorter times (higher redshift). 

For experiments attempting to measure the variation of the light speed 
at the present time, the derivative of c{t) (eq 29 with ri^/a^ << 1) will be 
more useful: 

= l_^n-n, = -^ + nA. (67) 

t=to ^ ^ 

Notice that this fraction is negative when matter dominates, and goes from 
zero at zero density to —1/2 at the critical universe density. A vacuum 
energy density opposes the gravitational effect of matter; when it dominates, 
the slope is an increasing function of time. 



1 dc 

c dt 



9 Underlying physics 

It really should not surprise us that gravity acting in concert with the cos- 
mological principle should cause a cosmic change in the light speed. It is well 
known that a rod stretching a short distance along the x direction in a grav- 
itational field g acquires a light speed of co(l -|- gx/c^) (see MTW[9]pl73). 

The effect of gravitational potential on light speed is also demonstrated 
by the Schwarzschild coordinates, where the coordinate light speed as well 
as the time on clocks are changed by the gravitational potential at a distance 
R from a central mass M. 

In Appendix C for FLRW we show that for a flat universe ($1 = 1) with 
the presently derived variable light speed, there is a gravitational field in 
the physicality region that increases linearly with distance from the origin. 
If we insert into eq 119 the mass of the universe inside the radius R, Mq = 
AttpqR^/S at time to, we obtain 

GqMq 

9 = (68) 

the Newtonian expression for gravitational field at a radius R inside a sphere 
of uniform density. This is another indication that the T,R coordinates 
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are obeying special relativity laws near the origin because an accelerated 
particle in the rest frame of SR has the Newtonian acceleration [11]. Note 
that g < indicates an inward pull on the galactic points towards the origin 
of the stationary axis, which we can interpret as the cause for the universe 
expansion to slow down (for A = 0). 

However, the FLRW light speed is a function of t not R. Clearly, in 
the FLRW universe, light speed and gravitational red shift are dominated 
by the homogeneous and isotropic symmetry assumed for the universe and 
not solely by the gravitational fields that dominate the Schwarzschild light 
speed. 

Thus, just as the assumption of homogeneity and isotroopy equires the 
universe to be either expanding or contracting, it seems to require the phys- 
ical light speed to depend on this rate of expansion or contraction. 



10 To observe c(t) 

Current measurements with atomic clocks ([6], [7]) have achieved an accu- 
racy that indicate the frequency of atomic spectra do not change with time. 
(Of course, when measured on a frame moving at a different velocity or 
in a gravitational field, frequency docs change). On an inertial frame, this 
means that the fine structure constant a/ and the Rydberg constant RooC 
(expressed as a frequency) do not change with c{t){vin). 
The fine structure constant a/ in SI units [13] is 



=2 



and the Rydberg constant is 



2 m-ec e me , . 

R c = af-—-= 3 . (70) 

•' A-kU eQ(47r^)3 

Note that the Attcq is often omitted in the fine structure constant since it is 
unity in Gaussian coordinates, but it is essential here if we are to consider 
a variable c{t) . 

For these to remain constant while keeping e, fi and mc? constant re- 
quires that 

eo(i)c(t) = — ] =k = eo(to)co, a constant. (71) 
IJ,o[t)c[t) 

The most straight forward way to observe c(t) is to find a way to di- 
rectly measure the light speed or the atomic spectra wavelengths with the 
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same precision and stability that we can now measure spectra frequency. A 
fractional change in speed or wavelength should be 6xl0~^^ in 100 sees or 
2x10^"'^^ in a year if c/cq = {t/toY^^. With this much sensitivity, however, 
an observation would have to separate out the effects on light speed of the 
gravitational forces of local masses like the earth, the moon, and the sun. 

c{t) might affect all of distant observations. For instance in the measure- 
ment of supernova Ia[16] it will affect the measurement of acceleration of 
galaxies. It's possible that the whole apparent acceleration may be explain- 
able by c{t) instead of dark energy (A). This would require the luminosity 
of the super novae la to decrease with increasing c. For a fiat universe the 
apparent distance is given by 

dL = ^(1 + ^)[1 - (1 + z)-i/3][L(c)]i/2 (72) 

where L is the fraction by which the luminosity is changed by c{t) /cq . (Note 
that if Hq is the reciprocal of the measured slope of d-L vs z for small z, then 
H'q = 3Ho/2 and coto^^o = 3/4. In general, H^^ = iJo [1 + /2 - J) a] ) • For the 
observations to be entirely explained by c{t) instead of dark energy would 
require the luminosity to vary as = (l + z)~^/^ for a flat universe without 
dark energy. Recently, this apparent acceleration has been confirmed using 
distant clusters as a standard candle [17]. For this to be explained by c{t) 
would require it to have the same c{t) dependence. 

Other astronomical observations that might be affected by c(t) are cos- 
mic background radiation, gravitational lensing, and dynamical estimates of 
galactic cluster masses. 

Unfortunately, the c{t) calculated herein solves neither the flatness nor 
the horizon problem without inflation: The flatness problem changes lit- 
tle because the Hubble ratio has the same dependence on the universe scale 
factor a{t). The horizon problem remains because c(t) enters both the trans- 
verse speed of light and the radial speed of matter. At the time of the release 
of the CBR photons, without inflation light could have traveled laterally only 
/q* c3^t/ra = 1/(1 -|- z)^/^ = .04 radians for 17 = 1,2; = 3000, and so could 
not have interacted with matter separated by more than this angle. This 
would allow only .002 of the universe to have interacted. 

11 Conclusions 

From the cosmological principle of spatial homogeneity and isotropy we can 
obtain the FLRW metric that describes a universe of inertial frames attached 
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to expanding galactic points with FLRW differential comoving coordinate 

times the scale factor a(t) interpreted as a physical differential distance. The 
FLRW metric is Minkowski-like in its radial derivative. Locally, SR applies, 
so a stationary rigid frame attached to the origin has a Minkowski metric. 
Thus, for a radial world line we can use a Lorentz transform from FLRW 
to the stationary frame that keeps the two Minkowski world line elements 
invariant in order to obtain time and distance coordinates to describe radial 
movement in the universe close to the origin. Because the FLRW metric has 
a time varying coefficient multiplying the space differential, this produces a 
velocity between the galactic points and the stationary frame that is a func- 
tion of time and distance. This places a double constraint on the stationary 
frame variables. A consistent limiting process to zero distance from the ori- 
gin with codt replaced by c{t)dt satisfies this double constraint by a variable 
light speed c{t) oc y^da/dt, the square root of the rate of change of the scale 
factor of the FLRW universe. By homogeneity, the origin can be placed 
on any galactic point, so that this variable light speed enters physical laws 
throughout the universe. The replacement of the FLRW differential time dt 
by dt = c{t)d.t retains the FLRW solution of the Einstein field equation, but 
with a variable light speed physics. 

Although three orthogonal rigid axes are inadequate to describe three 
dimensional motion in accelerating fields, it is possible to describe one di- 
mensional motion on a single axis. We have done this for the FLRW universe 
by finding radial covariant transforms from FLRW for all distances whose 
differentials remain close to SR Minkowski with this same variable light 
speed out to a red shift of 2 for a flat universe. 

I have shown that the physical coordinates on the stationary frame near 
the origin have a gravitational field for a flat universe that increases linearly 
with radius just like the Newtonian field for a spherical distribution of uni- 
form mass density. Like Schwarzschild, a gravitational red shift is predicted 
for a distant stationary light source observed at the origin of the FLRW 
universe. 

To summarize, I am persuaded that the physical light speed throughout 
the FLRW universe is proportional to y'da/dt because (1) in the limit of 
zero distance from the origin a radial Lorentz transform from FLRW to a 
stationary rigid frame requires it; (2) all radial covariant transforms from 
FLRW coordinates that I have investigated that are SR Minkowski close to 
the origin have this same variable light speed; (3) we can use the Einstein 
field equation to calculate the transformed distance vs time for galactic 
points and light that behave in a physically sensible way; (4) the transformed 
gravitational field in the physicality region for a flat universe is Newtonian 
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for a spherical distribution of uniform mass density and can be considered 
the cause of the deceleration of the universe (when dark energy can be 
neglected) . 

Just as the assumption of homogeneity requires the universe to be either 
expanding or contracting, it seems to require the light speed to depend on 
this rate of expansion or contraction under the influence of gravity. 

One of the radial covariant transforms from FLRW has a distance co- 
ordinate that remains physical for all distances. We can interpret this to 
be a global reference distance (used in Figs 1-5), although the time of this 
transform becomes unphysical at large distances. 

Some other physical "constants" that depend on the light speed must 
also be changing with cosmic time. I have suggested some constraints on 
this variability: (1) keeping constant the rest mass energy, the gravitational 
energy, and the Schwarzschild radius, and (2) keeping frequency of atomic 
spectra constant, which means the fine structure constant, and the Rydberg 
frequency. These still make possible the geometrization of relativity with an 
adaptation of vectors and tensors such as the energy-momentum vector, the 
stress-energy tensor, and the electromagnetic field tensor. 

This c(t) should be observable by direct measurement of light speed 
or spectral wavelength if they could be measured to the same precision as 
frequency, and the if effects on light speed of the gravitational forces of 
nearby masses like the earth, the moon, and the sun could be isolated. 
It should have an impact on understanding distant cosmic observations. 
Perhaps it will provide an alternative to dark energy to explain the apparent 
acceleration of galaxies via supernova la. Analysis of cosmic background 
radiation, gravitational Icnsing, and dynamical estimates of galactic cluster 
masses could also be affected. 

But the recognition of this c(t) does not solve the flatness nor horizon 
problems without inflation. 
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Appendix 

A Transformed diagonal coordinates with physi- 
cal time 

A.l Pcirtial differential equation for V — V/c{t) 

We will be considering radial covariant transforms for diagonal coordinates 
that eq 44 makes 

V = -^. (73) 
aTt ^ ' 

For diagonal coordinates with physical time at all t and ^ = 1- Thus, 

eq 42 becomes 

Tt = 7- (74) 

This automatically guarantees the Lorentz time dilation {dT/dt)^ = I/It = 
1/7 (eq 40). We need only find a transform for which B ^ 1 close to the 
origin. 

We proceed by finding a differential equation with V as the only de- 
pendent variable. Thus, we can write a formula for T, using eq 74 and eq 
73: 



a 



T = t+ T^dxt = t+ {—jV)dxt, (75) 
Jo Jo c 

where we have used the boundary condition that at x = 0, T = t, and the 

symbol dxt signifies integration with x a-t constant t. It can be partially 

differentiated with respect to t (giving 7) and then with respect to x 

with the use of eq 33, noting that d'y = j^VdV and 1 + V^j'^ = 7^, we 

obtain a PDE for V: 




A. 2 The solution for V, R. and T for all a without using 

c oc d 

Eq 76 can be rewritten as 



V{1 - f 2) 



(77) 



c 



where the subscript on the partial differential indicates the variable to be 
held constant. This can be integrated with an integration constant Iuk. 
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Since the integration is done at constant R, then k = k{R), and inversely, 
R = R{k). Integrating eq 77, we get 

V = -^^=, (78) 



\/5 



where the sign of k will be positive for an expanding universe, where the x 
points will stream out radially past a point at R. 

At this point, R is an unknown function of k. The various possible 
coordinate systems which solve our PDEs are characterized, in large part, 
by the function R{k). But for all, in order for V to vanish when R = (see 
definitions sect 1), k must also; so always 

k(0) = 0. (79) 

We note that as long as k{R) remains finite, V goes to —1, and V goes to 
—c{t), for a{t) = 0, i.e. for t = 0, the horizon. 

Let us now look at lines of constant k{R), i.e. constant R, in t,x space. 
Eq 33 can be integrated for x with use of eq 78 at constant k to give the 
following: 

For an expanding universe, we have set the upper limit at oo, because we 
expect that if R is kept constant the galactic point x that will be passing any 
given R will eventually approach zero as FLRW time t approaches infinity. 

At this point, we have obtained V = V*{t, k) from eq 78 and have also 
obtained the function xit,/^)- We can in principle invert eq 80 to obtain 
K in terms of t and x- = K{t,x)- This gives us the velocity function 
V{t,x) = V*{'t, K{t,x)- If the function R{k) were known, we would then 
also have R{t, x) = R{K{t, x))- 

T{t, x) can be found by noting from eqs 73 that 

aV^ aV 
T^ = -—Tt = -—j = K. (81) 

By substituting eq 81 into eq 75, and integrating over k instead of x by 
dividing the integrand by the partial of eq 80 with respect to k, we find an 
expression for T{t,x)' 



T{t,x) = t + 



1 



ds^, (82) 
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where k is put equal to K{t, x) after integration at constant k in order to 
get T{t,x). 

This completes the solution. Since k(R) can be any function that van- 
ishes at the origin, there thus exist an infinite number of solutions for our 
transformed coordinates with A = 1,C = 0. 

A. 3 Independent determination of c{t) 

To determine physicality, we will next find 1/Rx close to the origin. R^. = 
K^/ k'{R) can be written in an inverted form by taking the derivative of eq 
80 with respect to k at constant t: 



— = k' 



dK 



dy^Jt ^ a2(l + ^)3/2 



To be physical B = ^a/R^ ^ 1 as R approaches 0. Putting 7 = 1 and 

k(0) = in eq 83, multiplying by a, setting the resultant equal to one, and 
changing the integration variable from t to a{t) gives 

'=^'(0)/°°!^. (84) 

remembering that the dot indicates differentiation by t. k'(0) is a constant 
to be determined by c(to) = 1- Note that the integral of eq 84 is independent 
of the functional form of k{R), and is therefor the same for all k(R). It was 
Eq 83 that gave me the first indication that the fight speed (cR^/ja) was 
variable, and that it was the same near the origin for all k{R), which would 
be a requirement that it was indeed the physical light speed. 

Eq 84 is an integral equation for c(i) . By differentiation of both sides of 
eq 84 by a, we can obtain 

^(*) = (8^) 
k'{0) 

which, as we should expect, is the same c(t) (see eq 29) we showed for all 
physical coordinate systems for k'(0) = d{to)/co = ao-ffo/co- This indepen- 
dent derivation of c{t) confirms the validity of carrying the series expansion 
to second order since complete transforms give the same c{t). 

A. 4 Zero density universe = 

It is interesting to consider the limiting case of a zero density universe: 
fl = 0,flr = 1, aoHo = 1 (eq 57). Eq 29 makes c = 1. Eq 60 makes 
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a = Hq for all t,x- Integrating gives a = t. Eq 80 gives x = cscli~^t/«;, or 
K = K{t, x) = i sinh x- We can then find from eq 78 that V{t, x) = — tanhx 
and from eq 34 that 7 = cosh x so that 

c=l = ^ = - ^ = —. (86) 

7a CLK jt ClK 

Thus the physicality condition is met for all R with R = K and A = 1,B = 1, 
so that the complete transform with eq 82 becomes 

R = tsmhx, .„„x 
r = tcoshx. ^ ^ 

These coordinates have been known ever since Robertson [18] showed that 
this transformation from the FLRW co-moving coordinates at zero density 
obeyed the Minkowski metric. What is new is that this solution was derived 
from the equations we obtained for our physical time transforms with A = 1. 
It can also be obtained from the physical distance transforms (B = 1) since 
eqs 76 and 92 for V become identical with = and a/c = a = t. It is the 
only known rigid physical coordinate system for all times and distances in 
a homogeneous and isotropic universe. R is plotted vs T in Fig. 5 to show 
how similar it is to the physicality region of Figs. 2-4. 

B Transformed diagonal coordinates with physical 
distance 

B.l Pctrtial differential equation for V. 

For diagonal coordinates with physical dR for all t and X) -B = 1, so eq 43 
becomes 

R^ = a-f (88) 

By integration we find 

R = ar^dxt, (89) 







and partial differentiation with respect to t gives 

rx rx 
Rt = ca -fdxt + a 7t^Xt- (90) 
Jo Jo 

We can then find V from eqs 39, 88, and 89 as 



C7 C7 



X 

cd I -fdxt + a I ItOxt 



:d I -ydxt 
Jo 



(91) 
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This is an integral equation for V. It can be converted into a partial differ- 
ential equation by multiplying both sides by 7 and partial differentiating by 

V.+ '^VVt\=-a = -'y-^. (92) 
c J c at 



1' 



Note that this is substantially different from the eq 76 for V that we ob- 
tained for physical time This means that it is not possible to find diagonal 
transforms with both physical time and physical distance for all values of t 
and X (except for = 0). It is possible to have either one or the other be 
physical at all t and x with the other being physical only close to the origin. 

B.2 General solution for V 

Eq 92 can be solved as a standard initial-value problem. Let W = —V. Eq 
92 becomes 

a 1 da , n, , , 

W^- -WWt= - — {l-W^) (93) 

Define a characteristic for W{t,x) by 
and 

(The subscript here indicates differentiation along the characteristic) . If we 
divide eq 95 by eq 94 we get 

,d_W _ Ida jl-W^) 
^ dt^"- adt W 

This can be rearranged to give 

W{dW)c {da)c 



- 1 



(97) 



This can be integrated along the characteristic with the boundary condition 
at X = that W = Q and a = Qc- 



a" 



l-W^ = ^ (98) 
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This value for W can be inserted into eq 94 to give 



:|^)c = A/i-"J (99) 



We can convert this to a differential equation for a by noting that cdtc 
die = Tidac 



a? 



or for a using eq 60: 



(^)c = -oaWl- ^> (100) 



{-^)c = -aoHoa^E{a)Jl-^^ (101) 

This can be integrated along the characteristic with constant Oc, starting 
with a = Oc at X = 0. This will give x = ^if^, Q^c)- This can be inverted to 

obtain ac{c(, x). When this is inserted into eq 98, we have a solution to eq 
93 for W{a,x)- Then W{t,x) can be found from the relation 

r da 

t= — (102) 
Jo a 

So far, I have not used the relation c = coaE{a). To get W{t,x) we need 
to use it: 

r da 1 r da , , 

t= — = —rr I -^r^ (103 
Jq ca cqHo Jo a^E^ 

B.3 Obtaining T, R from V 

Eqs 33, 39, and 73 show that 

»' = -5(|)„ = 2|t = fa (104) 

c at c aTt 

Thus T has the same characteristic as W so that {dT/dx)c = and T 
is constant along this characteristic: 

T{t, x) = T{te, 0) = tc = t{ae{t, x)) (105) 

where t{a) is given in eq 103. This gives us the solution for T{t,x)- 
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The solution for R can be obtained by integrating eq 89, using 7 from 
eq 98 and ac{t,x) from eq 101: 

rx rx 

R{t,x) = a 7dr]t= ac{t,r])dr]t (106) 
Jo Jo 

Alternatively, for ease of numerical integration we would like to integrate 
dR along the same characteristic as T and W. This can be obtained from 
the PDE 

{dR/dx)c = Rx + Rt{dt/dx)c (107) 

If we insert the values for these three quantities from eqs 88, 104, and 94 , 
we get 

{dR/dx)c = 70 + {cW/a)^a{-aW/c) = a/7 = a^/ac- (108) 



C Similarity solutions for flat universe, = 1 

I have found similarity integrations for the special case of f2 = 1 where 

o = aQit/tQ)^/"^ and c = (tg/i)^^^- To simplify notation let us normalize 
time io t/tQ t, a/ aQ ^ a, and xoo/coio x, T/to — > T, R/cqIq — > R, and 
let W = -V. 



C.l Physical Distance 

Eq 92 then becomes 

W:c - t^'^WWt = ^t-^'^{l - W^) (109) 

This can be converted into an ordinary differential equation (ODE) by let- 
ting 

u = x/t^/'^ (110) 

so that eq 109 becomes 

W'{l + '^) = l{l-W% (111) 

where the prime denotes differentiation by u. 

Similarly we can find ODE's for T and R by defining: 

T/t = q{u), (112) 

and 

R/t^/'^ = s{u), (113) 
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where q{u) and s{u), from eqs 73 and 39, are given by the coupled ODE's: 

q\l+'^) = qW, (114) 

and 

AW+l) = ls (115) 

It is useful to find that q = , s' = ^, s = ^{u + 4:W)/3, and A = j/Tt = 
(1 + uW/A) /r, so T = jH, and R = i3/4^(n + 4W) /3. 

For small values oi u, W = u/2, q = 1 + s = u, and R = t^f^x = 

ax. We also note that — 1 + confirming that these coordinates 

have physical time close to the origin, justifying c(t) = t~^f^. An alternate 
approach would be to start with c{t) unknown, but of the form c = t~^. This 
will also yield a similarity integration of eq 92 with Tt = l + 2bW'^ for small 
u. In order to make the time physical close to the origin, Tt has to equal 
1 + so that b = 1/4, confirming that for these diagonal coordinates to 

be physical at the origin c oc da/dt. 

As t ^ 0, n ^ oo, and ^ 1 - 760/(« + 4)^ q ^ {u + 4)'^/760, 
and s — ^ .0121(tt + 4)^. T and R both remain finite at this limit with 
T .00132x4, and R .0121x3. r^-^^ f^^^ ^^^^^ ^^^g 

zero when 

t goes to zero results from synchronizing T with t sX t = 1 and not at i = 0. 

The distance R and time T can be found from the numerical integration 
of the coupled ODE's. The paths of galactic points are those for constant x. 

The path photons have taken reaching the origin at ti is found by calculating 
Xp vs t and using the transform to T, R. Thus, for = 1 

rt 



Xri 



[ -dt=A{ty'-t^/^) (116) 

•7*1 



For light arriving now, ti = 1, the value of Up becomes 

«P = 4(^-l) = 4(c-l), (117) 

where we inserted c = to obtain the relation of c to Up. 

Galactic and photon paths are shown in Fig. 1 and 2. An approximate 
upper limit of physicality is shown by the heavy dotted line: A = .953, 
= 0.253, u = 2.0, R = 2.30^-^1 = 1.35r3/4. R vs T at t = provides 

a non-physical horizon: R^ = 1.747T^^^. 

It is also interesting to calculate the acceleration g. If we insert the 
values of V, R, and a/c in eq 142, we obtain 

1 



-9 



u 



87t5/4 [l + uW/4. 



(118) 
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where the units of g are co/to- For small u close to t = to, g goes to zero as 
—u/8. Since this is the region with physical coordinates, it is interesting to 
express this in unnormalized coordinates: 



-9 



Icq R 
8 to Coto 



GqPo-^R, 



(119) 



where we have used $7 = 1 in eq 56. For small t, —g goes to oo as 1/ {2jt^^^) = 
1.2t~^^'^ along the light path. At the physicality limit, —g = .69co/to = 
12.5 X lO-^m/sec^ 

g can be obtained from a gravitational potential using g = —d(f)/dR, 
which for close distances is: 



= Cn 







'16" 



(120) 



The slope of the light path in Fig. 1, a coordinate velocity of light, can 
be shown in normalized units for this incoming light path to be 



,dR, 



-(1 + u/4) 



l-W 
1 + W 



.1/2 



For small u , 



(121) 



VL^ - (1 - U/^)- 



For the outgoing light path 



^l = (^)l = (1-«/4) 



l + W 
l-W 



1/2 



For small u. 



(122) 



(123) 



(124) 



Thus, the coordinate light speed has a different u dependence and therefore 

R dependence for the two directions of the light path. This differs from the 
Schwarzschild solution that has the same R dependence of the coordinate 
light speed for both directions of the light path. 

The observed light at the origin v that is emitted from a stationary 
source at R as uq is also smaller than the same light emitted at the origin. 
This can be shown to be 



- = {-)n=- 
uo ^ dt jA 



q 



1 - 



(125) 
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Close to the origin it is: 



V 



1 



y2 

2 



1 



u 



1 



2(p 

'72 ' 



(126) 



This, of course, is the same as a dilation effect for a co-located galactic 
point at R that shows up as a gravitation red shift at the origin due to 
homogeneity of t. This is interpreted in Sect 9. 

C.2 Physical Time 

There is also a similarity solution for physical time, A = 1, for Q = 1. With 
the same normalizations as above, using eq 76, the ODE for W is 



W'{W + 



-W{l-W^), 



(127) 



with the ODE's for T, R, q, s, Xp) the same. This is the same as physical 
distance for small u, but differs numerically at large u. Useful relations for 
physical time are obtained from the general solution in Appendix A. R = 2k, 
q = 7(1 + uW/4), s = 2-fW, and B = -f{u + 4W)/3s = 2(1 + u/AW)/3. 
These can be used to find the gravitational field from eq 142: 



9 



2*5/4 u + AW 

and the coordinate velocity for an incoming light path: 



VL = -; 



1 + u/4W 



1-W 
1 + W 



1/2 



(128) 



(129) 



Eqs 128 and 129 approach the same values as physical distance for small u. 



D Gravitational field in the transformed system 

Wc wish to find the components of the radial acceleration of a test particle 
located at R in the covariantly transformed system. We will do this by 
calculating the FLRW components of the acceleration vector and find the 
transformed components by using the known diagonal transforms. For the 
FLRW components, we will use the metric 

ds'^ = - a^dx^ - a^r'^dd^ - a\'^ sin^ edcj)'^ . (130) 
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Let 

= X, = f, =4>, x^ = t, (131) 
and the corresponding metric coefficients become 

344 = 1, 511 = -a^, 922 = -a^r^, 933 = -a^r^sin^f. (132) 
For any metric, the acceleration vector for a test particle is 

A^=^+ T^^U^U", (133) 

where the F's are the affine connections and is the velocity vector of 
the test particle. In our case the test particle is at the point R on the 
transformed coordinate, but not attached to the frame so that it can acquire 
an acceleration. Instantaneously, it will have the same velocity as the point 
on the transformed coordinate, and its velocity vector will therefor transform 
the same as the point (eq 38). 

We will be considering accelerations only in the radial direction so that 
wc need find affine connections only for indices 1,4. The only non-zero 
partial derivative with these indices is 

Is = (m 

The general expression for an affine connection for a diagonal metric is 



dx^ ^x^^ dx^ 
The only three non-zero affine connections with 1,4 indices are 



(135) 



rfi = ad, rii = ri4 = -. (i36) 

The acceleration vector in FLRW coordinates of our test particle moving at 
the same velocity as a point on the transformed frame becomes 



A>^ = ^+rii(C/4i7i + C/ii7^), 
Using and in eq 37 we find 



(137) 



(138) 
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Since the acceleration vector of the test particle at R in the transformed 
coordinates will be orthogonal to the velocity vector, it becomes 



Ds 



(139) 



is the acceleration of a point on the R axis, so the gravitational field 
affecting objects like the galactic points is the negative of this, g is defined 
so that mg is the force acting on an object whose mass is m (eq ??). For 
a range of time in which c(t) is reasonably constant, g = the normal 
acceleration. Since the vector A'^ will transform like dT,dR (eq 7): 



iR 



-RtA^ + RyA^ 
c 



so that 



7^ 



dt a 
/ R 



^Rt + 
c 




With the use of eq 39, this can be simplified to 



R^ 



dV 



'-V 



(140) 

R^ 
(141) 

(142) 



In terms of the normalized coordinates for a fiat universe (Appendix C), 
this becomes 



9 



s 

'JJl 



w 



(143) 



The acceleration g can be thought of as the gravitational field caused by 
the mass of the surrounding galactic points, which balances to zero at the 
origin, where the frame is inertial, but goes to infinity at the horizon. It is 
the field which is slowing down the galactic points (for A = 0) . It is also the 
field that can be thought of as causing the gravitational red shift (App. C). 
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Figure 1: Physical distances (i?/coto) for = 1 plotted against the normal- 
ized time on clocks at the origin (t/^o) for various galaxy paths (labeled by 
their red shift z) and for the light path which the photons take after emission 
by any galaxy that arrives at the origin at t/to = 1- Notice that the slope 
of this light path close to the origin is cq = 1, where the coordinates R and 
t are both physical. The light path starts at the far horizon at t = 0, not at 
i? = 0, but the transformed time T on the stationary frame remains finite 
(see Fig. 2). The dashed line in Fig. 1 shows the approximate upper limit 
of physicality, where both R and the transformed time T are physical. 
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Figure 2: Physical distances (ii/coto) for = 1 vs the transformed time 
(T/to) on clocks attached at i? for various galaxy paths (labeled by their 
red shift z) and for the light path that arrives at the origin at T = to = 
.75{1/coHq). The horizon is the locus of points where t = 0. The heavy 
dotted line in Fig. 1 shows the approximate upper limit of physicality for 
the transformed coordinates {A = .95, that is < 5% error in physical time 
rate Tt). The slope of the light path is very close to c{t) out to the limit 
of physicality. Light is emitted at finite T allowing transformed time for 
galactic points to move out from R = before emitting their light. 
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Figure 3: Physical distancc(i2/coto) lower density universe (17 = 
1/2, 17^ = 1/2) plotted against the transformed time (T/to) on clocks at- 
tached at R for various galaxy paths (labeled by their red shift z) and for 
the light path that arrives at the origin at T = to = •767(l/co-ffo)- '^^^ 
horizon (t = 0) and the physicality line i^A = .96) occur at later times and 
shorter distances than for a flat universe (Fig 2), but not as much as for the 
empty universe (Fig 5). Similarly, the light path is straighter than Fig 2, 
but not as straight as Fig 5. 
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Figure 4: Physical distance (i?/coto) vs. transformed time (T/coto) foi' dark 
energy = 1/4, Oa = 3/4). The present time to = A07{1/coH'q). Notice 
the inflection points on ah curves where the dark energy density becomes 
larger than the matter density. 
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Figure 5: Physical distancc(i?/coto) for the empty expanding universe {Tl = 
0, Or = 1) plotted against the transformed time {T/to) on clocks attached 
at R for various galaxy paths (labeled by their red shift z) and for the light 
path that arrives at the origin at T = t^. The horizon is the locus of points 
where t = 0. All lines are straight and physical, since there is no space 
curvature, and the light speed is c{t) = cq. The remotest galactic point 
travels from the origin at T = out to coto/2 at the light speed cq. 
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